Coeovolutionary Threshold Dynamics 
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We present a generic threshold model for the co-evolution of the structure of a network and the 
state of its nodes. We focus on regular directed networks and derive equations for the evolution of the 
system toward its absorbing state. It is shown that the system displays a transition from a connected 
phase to a fragmented phase that depends on its initial configuration. Computer simulations are 
performed and confirm the theoretical predictions. 
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The relation between a cause and its effect is usually 
abrupt in complex systems, in the sense that a small 
change in the neighborhood of a subsystem may (or not) 
trigger its reaction. This mechanism is at the heart 
of many models of self-organized criticality [1] where a 
cascade starts when the system has been frustrated be- 
yond some threshold, e.g. the angle of a sand pile, but 
also in models for the diffusion of ideas in social net- 
works [H, 0, 0, S] where the adoption of a new idea 
requires simultaneous exposure to multiple active ac- 
quaintances, and in integrate- and- fire neuron dynam- 
ics [i[ where the voltage on a single neuron increases 
until a specified threshold is reached and it suddenly 
fires by emitting an action potential, thereby quickly 
returning to its reference. These types of model con- 
sist in cascading propagations on a fixed topology, i.e., 
a network of some sort, until a frozen configuration is 
reached, but they do not incorporate the well-known 
feed-back existing between network topology and dynam- 
ics 0, H !, ED, El EH, El, El El , namely that the topol- 
ogy itself may reorganize when it is not compatible with 
the state of the nodes. This reorganization may originate 
from homophily and social balance in social networks or 
synaptic plasticity in neuron dynamics. 

The purpose of this paper is to introduce a model for 
coevolutionary threshold dynamics (CTD). Let us de- 
scribe its ingredients in terms of diffusion of opinions in 
social networks [l6|, E3 while keeping in mind that the 
model is applicable to more general systems. The system 
is made of a social network of interaction, whose N nodes 
are endowed with a binary opinion s, + or — . The dy- 
namics is driven by the threshold 0, such that < <\> < 1 
and, in most cases of interest, > 1/2. At each step, 
a randomly selected node i evaluates the opinion of its 
hi neighbours. Let fa be the fraction of neighbours dis- 
agreeing with i. If <pi < (j), node i breaks the links toward 
those disagreeing neighbours and rewires them to ran- 
domly selected nodes. If fa > 0, i adopts the state of the 
majority. By construction, the dynamics perdures until 
consensus, i.e., all agents having the same opinion, has 
been attainted in the whole system or in disconnected 
components. This absorbing state obviously depends on 
the threshold <p but also, as we will discuss below, on its 
initial condition. 



fraction of disagreeing neighbours 



FIG. 1: Update process of CTD for two different configura- 
tions of neighbors. When one out of four neighbors is in a 
different state, the central node breaks its links and creates a 
new link to a randomly chosen node. When three out of four 
neighbors are in a different state, the threshold <j> is exceeded 
and the central node thus adopts the majority state. 



A complete analysis of CTD requires extensive com- 
puter simulations, which is not the objective of this pa- 
per. We will instead focus on a simplified version of the 
model that can be studied analytically and pinpoint the 
key mechanisms responsible for its behaviour [27]. In 
this simplified version, the network is directed and all 
the nodes have two incoming links, i.e. each is influenced 
by exactly two nodes, while their out-degree is initially 
Poisson distributed. Moreover, we will take <\> = 1/2, 
such that CTD now simplifies as follows. At each time 
step, a node i is selected at random. If i is surrounded by 
two nodes with different opinions, it switches its opinion, 
i.e., Si — > — Si. If the opinion of only one of its neighbours, 
say j is different, i cuts its link from j and reconnects to 
a randomly chosen node, i.e., its in-degree remains con- 
stant. It is interesting to stress that the choice 0=1/2 
for a directed network with a constant in-degree of two 
corresponds to the unanimity rule [18| when no rewiring 
is implemented. This model is well-known to exhibit a 
non-trivial relation between initial and final densities of 
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+ nodes, denoted by n +; o and n +;oo respectively. We 
will show that the addition of the rewiring mechanism 
leads to a transition from a connected phase with con- 
sensus where all the nodes asymptotically belong to the 
same cluster, to a fragmented phase where two discon- 
nected clusters of different opinions survive. The critical 
parameter of this transition is shown to be the initial 
density n+ ; o of + nodes, i.e., 

n +; o < n c or n +; o > 1 — n c one single component, 
n c < n+ ; o < 1 — n c two disconnected components, (1) 

where n c is the critical density. 



In order to analyze the system dynamics, let us follow 
the approach proposed in [18] and focus on the number 
N SQ . SlS2 of configurations where a node in state so re- 
ceives its incoming links from a node in state si and an- 
other node in state S2- Let us denote by {so; such a 
triplet of nodes. By construction, si may be +1 or —1 and 
^s s 1 s 2 Ns ;s 1 s 2 ~ N . Moreover, the order of the links 
is not important and therefore N SQ . SlS2 = N SQ . S2Sl . By 
neglecting higher order correlations than those included 
in N So;SlS2 , it is possible to derive the set of equations 
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where n + and n_ are the density of + and — links respec- 
tively. 7r^ >_ (tt , + ) is the probability for a randomly 

selected + (— ) node to switch its opinion to - (+). By 
construction, this quantity is the probability that a ran- 
dom + (— ) node is connected to two — (+) nodes 
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where N+ = 
total number of + and — nodes respectively. 

Let us describe in detail the first equation for 7V +;++ , 
the other ones being obtained in a similar way. Its evo- 
lution is made of several contributions. The first term 
is the probability that a { — ; ++} triplet is selected and 
transforms into {+;++} by unanimity rule. The second 
term is the probability that a {+; H — } triplet is selected 
and the rewired link (originally from + to — ) arrives on 
a + node (with probability n + ). The last two terms ac- 
count for possible change of the state of one of the two 
neighbours in the triplet, as they may also switch their 
opinion because of a unanimity rule in another triplet, 

and are evaluated by using the aforementioned 7r^ >_ 

and 7r >_(_. 

As discussed in [l8|, several initial conditions may in 
principle be chosen for the system of equations ([2]), each 
of them leading to its own trajectory in the 6-dimensional 
dynamical space. Such initial conditions are subject to 
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are the total number of + (— ) incoming neighbors in 
the triplets. Relations (jlj) simply mean that each node 
i that is a neighbor in a triplet {s x ; SiS y } is also at the 
summit of another triplet {s^; s x t s y > } (as it also receives 
two incoming links by construction). In order to select 
one of the several configurations N So;SlS2 that still satisfy 
the above constraints, we will further assume that the 
initial configuration is uncorrelated and therefore that 
each node has the same probability n +; o to be +. Among 
all the possible configurations for which A+ = An +; o, we 
therefore select the initial condition 
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Simulations 
Numerical integration of Eqs. 2 
Order parameter r 
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FIG. 2: In the upper figure, we plot the relation JSJ) between 
the initial density and the final density of + nodes, evaluated 
by integrating the set of Eqs.((2]) and by performing numerical 
simulations of a network made of 1000 nodes. We also plot 
the order parameter r=(l/2 — |l/2 — n+,oo|) which confirms 
that the system actually breaks into disconnected components 
when n c < n+,o < 1 — n c , i.e., we only plot © in the interval 
[0, 0.5] as the curves are symmetric around the point (0.5, 0.5). 
In the lower figure, we plot the probability density p(n+ ;0 o) 
that the absorbing state has a density n+ ;oo of + nodes for 
N 100 and N 1000 when n +;0 0.4. 



Before going further, it is instructive to look at the 
total number iV+ = ^2 Sl S2 N+ ;SlS2 of + nodes whose 
time evolution is obtained by summing over the first three 
equations of ([2]) 
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This relation shows that 7V_. ++ = 7V + .__ at stationarity. 
A careful look at the second equation of ([2]) shows, how- 
ever, that 7V +; __ has to decay until it reaches zero. The 
third equation of ([2]) also shows that the only stationary 

solution of 7V_|_ ; _| is also zero when iV_ ;++ = iV +; __ = 0, 

thereby confirming that the dynamics asymptotically 
reaches a frozen configuration where consensus is reached 
among connected nodes. The dynamics is therefore 

driven by two types of triplets: the triplets {+; } and 

{ — ;++} drive the system toward consensus, while the 
configurations {+; H — } and { — ;H — } allow for a topo- 
logical rearrangement of the network. This rearrange- 





FIG. 3: Visualisation of the initial and final states of one real- 
ization of the dynamics for a network made of N = 100 nodes. 
The initial density of + nodes n+,o = 0.4. The asymptotic 
network is made of two clusters. The + cluster is now made 
of 35% of the nodes. 



ment implies that the only frozen states are those cor- 
responding to consensus (in one or several clusters) [28[ 
and drives the division of the system into disconnected 
clusters. The competition between these two types of 
mechanisms is crucial for the transition (pQ). One should 
also note that models for opinion dynamics are known 
to exhibit coexistence of different opinions when applied 
to a static underlying network with modular structure 
|20|| . In the case of CTD, in contrast, it is the rewiring of 
the links that reorganizes the system into modules and 
thereby allows for coexistence. 

By integrating recursively the system of Eqs. ([2]) start- 
ing from the initial conditions ([6]), we obtain a non-trivial 
relation 



™+,oo(™+,o) 
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between the initial density and the final density of + 
nodes. This numerical integration confirms the above 
discussions, and clearly shows that a transition occurs 
at n c ~ 0.22 (see Fig. 2). One should insist on the fact 
that this relation differs from the standard exit probabil- 
ity measured when the dynamics takes place on a static 
network [2l|, [22|, [23[. In the latter models, n +)00 (n +) o) 
would measure the probability to end in a + consensus 
(in the whole system) starting from some initial density 
of + nodes. Relation J5J) is more reminiscent of the stan- 
dard unanimity rule [18| without rewiring, where the sys- 
tem asymptotically reaches a frozen state different from 
consensus at each realization, and where n +500 is the av- 
erage number of + nodes in this frozen state. Because 
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of the rewiring process, however, a non- vanishing value 
of n +>00 also implies that the system has split into two 
disconnected clusters and that a different consensus has 
been reached in each cluster. 

We have verified the accuracy of our calculations by 
performing numerical simulations of the model (see Fig. 
2). To do so, we have considered systems made of 1000 
nodes and have averaged the asymptotic density of + 
nodes (evaluated when the dynamics is frozen) over 1000 
realizations for each value of n+,o- In order to check 
that the system actually breaks into two clusters when 
n c < n +5 o < 1 — n c (see Fig. 3), we have also measured 
r = (1/2 — 1 1/2 — n+ j00 |). This order parameter would 
vanish if, for each realization, n+ 500 is either zero or one, 
while r = (n +500 ) if the system breaks into two clus- 
ters. The simulation results show an excellent agreement 
with the theoretical predictions and confirm the fragmen- 
tation of the network at the critical value n c . Finally, 
we have also looked at the probability density p(n +)OC ), 
i.e., Jq 1 p(x)dx = 1, that the absorbing state is made of 
n +500 7V nodes. This quantity is measured by performing 
5 x 10 4 simulations starting from the same initial condi- 
tion n +; o = 0.4. The distribution is shown to be peaked 
around its average, in contrast with the two delta peaks 
at and 1 that would be expected if full consensus had 
been the only absorbing state. 

To conclude, we have focused on a model for revolu- 
tionary threshold dynamics where the binary state of a 
node and its links coevolve. We have shown that the sys- 



tem may undergo fragmentation, a feature that has been 
observed in other coevolution network models, based on 
the Axelrod [l2[ or the Voter model [Hj] for instance, but 
also in the case of coupled maps with variable coupling 
strength [2J]. In the case of CTD, the critical parameter 
is the initial condition, as a sufficient fraction of + nodes 
is necessary for such nodes to survive and to separate 
from the main cluster. In this paper, we have focused on 
a simplified version of CTD where the underlying net- 
work is directed and regular, and where the in-degree 
has the smallest non-trivial value, i.e., two. Additional 
computer simulations are therefore required in order to 
explore the role of the threshold qb on the asymptotic 
state in more complex directed or undirected networks. 
Finally, let us point to an interesting generalization of 
CTD that would include two different threshold <j) r and 
(j) a for either rewiring links from disagreeing neighbours 
or adopting the state of the majority. Such a model 
would unify two seminal threshold models, namely the 
Granovetter model for the diffusion of cultural traits Q 
and the Schelling model for social segregation [25^ . 
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